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c/j , Dynamics of long-wave isocurvature perturbations during an in- 

flationary stage in multiple (multi-component) inflationary models is 
calculated analytically for the case where scalar fields producing this 
stage interact between themselves through gravity only. This enables 
to determine correct amplitudes of such perturbations produced by 
vacuum quantum fluctuations of the scalar fields during the multiple 
inflationary stage. Exact matching to a post-inflationary evolution 
that gives the amplitude of isocurvature perturbations in the cold 
dark matter model with radiation is performed in the case where a 
massive inflaton field remains uncoupled from usual matter up to the 
present time. For this model, isocurvature perturbations are smaller 
than adiabatic ones in the region of the break in the perturbation 
spectrum which arises due to a transition between the two phases of 
inflation, but they may be much bigger and have a maximum at much 
shorter scales. The case of an inflaton with a quartic coupling which 
remains uncoupled after inflation is considered, too. 



Abstract 







1 Introduction 



Inflationary cosmological models in which a de Sitter (inflationary) stage is 
produced by a number of effective scalar fields (inflatons) are called multiple 
(or multi-component) Jl| (see for a general review). A double inflation- 
ary model with two scalar fields ||-[f| is a specific case of them. Note that 
extended inflationary models or inflationary models in the Brans-Dicke 
theory of gravity may also be considered as belonging to this class of mod- 
els after transformation to the Einstein frame. Double inflationary models 
producing a step-like spectrum of initial adiabatic perturbations give a pos- 
sibility to reconcile the CDM model with observations without introducing 
neutrinos [L^]. If iV is the number of light scalar fields at the inflationary 
stage in such a theory (\m 2 \ -C H 2 , H = a/a, where a(t) is the scale factor 
of the FRW isotropic cosmologiocal model), then N independent branches of 
non-decaying quantum fluctuations of the scalar fields are generated during 
the inflationary stage similar, and in addition, to quantum fluctuations of 
gravitons (the resulting energy spectrum of the latter was first correctly cal- 
culated in [|T3J). However, only one linear combination of these fluctuations 
produces the growing scalar (adiabatic) mode which is usually assumed to 
be responsible for the formation of galaxies, stars (and other compact ob- 
jects) and the large-scale structure of the Universe. The other N — 1 modes 
are isocurvature fluctuations during the inflationary stage (they were first 
considered in fl4H ). 

Isocurvature fluctuations are less universal than adiabatic ones. First, 
they only appear in multiple, not single inflationary models. Second, they 
might not survive up to the present time (and typically do not). Really, they 
can exist now only if at least one of the inflaton scalar fields remains non- 
thermalized and uncoupled from the usual matter (radiation, baryons and 
leptons) during the whole evolution of the Universe from the inflationary 
era until the present period (so that the corresponding particles or prod- 
ucts of their decay constitute a part of cold dark matter) - a rather strong 
assumption. Finally, even so, their amplitude (in sharp contrast with adia- 
batic perturbations) does depend on the form of the transition from inflation 
to the radiation-dominated FRW stage. Nevertheless, isocurvature pertur- 
bations represent an interesting and important object of investigation, espe- 
cially because some candidates for such inflatons that might survive from the 
inflationary era up to the present time are already known - dilatons in the 
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Brans-Dicke theory and superstring induced theories (see |T3] for the latter), 



axions in "natural" inflation ||16|| , etc. 

So, in the present paper we consider isocurvature perturbations in the 
simplest case when N scalar fields have arbitrary self-interaction potentials 
but interact mutually through gravity only, i.e. the interaction potential 
V((fii, ...4>n) = J2\ V n ((j) n ). The general quantitatively correct expression for 
adiabatic perturbations generated in this model was obtained in First, 
we find the general solution for the behaviour of long-wave isocurvature per- 
turbations during the multiple inflationary stage and then determine the 
correct coefficients in it by exact matching to vacuum quantum fluctuations 
of the scalar fields in the approximately de Sitter background during the in- 
flationary stage (Sec. 2). After that, in Sec. 3, we investigate some cases 
where it is possible to match this solution to the radiation-dominated FRW 
model with a small admixture of non-thermalized massive particles ("cold 
dark matter"). The most interesting case with respect to cosmological ap- 
plications turns out to be the double inflationary model with two massive 
inflatons, the heavier one remaining uncoupled from usual matter after infla- 
tion. The isocurvature perturbation spectrum in this model has a maximum 
on small scales whose value may be rather large. 



2 Behaviour of perturbations during a mul- 
tiple inflationary stage 

We consider the following Lagrangian density describing gravity plus N scalar 
fields 

p N i 

L = -^ + ^(r 2 Mt-vA^)) (i) 

where fi = 0, ..,3, c = h = I and the Landau-Lifshitz sign conventions 
are used. Note that the N scalar fields interact only gravitationally. The 
space-time metric has the form 

ds 2 = dt 2 - a 2 (t)5 mn dx m dx n , m,n = l, 2, 3. (2) 

Spatial curvature may always be neglected because it becomes vanishingly 
small after the first few e-folds of inflation. The homogeneous background is 
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treated classically, it is determined by the scale factor a(t) and the N scalar 
fields 4>j{t). Their equation of motion is given by 

fl» = £!^(| + y jW ,)), (3) 

h + 3H& + V'(4>i) = 0, j=l,-,N (4) 

where a dot denotes a derivative with respect to t while a prime stands for a 
derivative with respect to <pj. From (f|), we get the useful equation 

N 

H = -AarCEfi. (5) 

i=i 

In these models, therefore, H always decreases with time. 

Let us turn now to the inhomogeneous perturbations. We consider a 
perturbed FRW background whose metric, in the longitudinal gauge, is given 
by 

ds 2 = (1 + 2§)dt 2 - a 2 (t)(l - 2ty)5 mn dx m dx n . (6) 

We get from the perturbed Einstein equations (exp(zkr) spatial dependence is 
assumed and the Fourier transform convention is $(k) = (27r)~ 3 / 2 J $(r)e~ lkr <i 3 k) 

$ = (7) 

AT 

$ + = AnG^faWj, (8) 

% + 3^4 + (^ + ^)^ = 40^-2^'$, j = l,..,N. (9) 

We see that when we have more than one scalar field, the dynamics of the 
perturbed system cannot be described by just one equation for the master 
quantity $. It is remarkable that without solving the system ([7|-[]) we can 
immediately write its two exact solutions desribing the adiabatic modes in 
the formal limit k — > 0: 

$ = Cjl-— fad?) + C 2 — , (10) 
\ a Jo / a 

5 -p- = ifd f ad? -C 2 ), j = l,..,N (11) 



-j a v jo 
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where a(t), </>j(t) satisfy the exact background equations (3,4) and G\ and 
C2 may still depend on k. The term with C\ is the growing adiabatic mode, 
the term with C2 is the decaying adiabatic one. The existence of these exact 
solutions directly follows from the observation made in |Tj| (see also the 
detailed explanation in 0) that there always exists a solution for scalar 
perturbations in the flat (JC = 0) FRW Universe which has the following 
asymptotic behaviour in the synchronous gauge in the limit k — > in terms 
of the Lifshitz variables: /i(k) = 3/i(k), A(k) = 0, 5<pj(k) = (with no 
dependence on t) irrespective of the structure and the properties of the energy- 
momentum tensor of matter. Knowledge of the solutions (10,11) is not, 
however, sufficient to find the amplitude of generated perturbations if the 
number of scalar fields iV > 1, in that case we have to integrate the system 
(7-9) completely in the limit k — > at the inflationary stage. 

Let us now consider a multiple inflationary stage with N background 
scalar fields being in the slow-rolling regime (N < N), N may depend on 
time. The energy density of all other scalar fields not being in the slow-rolling 
regime decreases exponentially with time and soon becomes negligible, thus, 
these fields should be simply omitted from the background equation (3). 
Then (3,4) simplify to 

H 2 = E^Wi), (12) 
i=i 6 

3#0, + H0;) = 0, j = l,..,N. (13) 

Now, for k <C aH, the system (7-9) can be solved in a way completely 
analogous to fl], ||. First, its solutions corresponding to growing adiabatic 
and non- decreasing isocurvature modes weakly depends on time, so for them 
Eqs. (8,9) take the form: 

* = yE^i- (14) 

SHSfa + Vj'Sfa = -2Vfi, j = l,..,N. (15) 
The general solution is 
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2H&g--d i ), i,j=l,..,N. (17) 



Here C\ and dj are integration constants, only N — l out of the N coefficients 
dj are linearly independent, and we will further use this freedom to add a 
constant term to them. The background quantities H(t),<f)j(t) are exact 
solutions of Eqs. (12,13). The mode with the coefficient C\ is the growing 
adiabatic mode as can be seen from the comparison with (10,11), the other 
N—l modes are the non-decreasing isocurvature modes. 

The expression for the decaying adiabatic mode immediately follows from 
the general expressions (10,11) which take the following form at the inflation- 
ary stage: 

• = c ^, ^ = _« (18) 

a <pj a 

for all scalar fields (including those which are not in the slow-rolling regime). 
The expression for decaying isocurvature modes of slowly-rolling scalar fields 
may be found, similarly to |J, by assuming that all quantities in Eqs. (7-9) 
are proportional to a~ 3 (t) multiplied by slowly varying functions of t (note 
that the approximate form (14,15) of these equations cannot be used now). 
The answer is 

$ = * = 0, 6cf> j = T ^— 53^ = 0, j = l,..,N (19) 

which may be easily verified by direct substitution using (12,13). Finally, 
all other 2(iV — N) scalar modes connected with non-slowly- rolling scalar 
fields are decreasing isocurvature ones, too. We shall consider them below in 
connection with matching to a post-inflationary era. 

Let us return to the most interesting non- decreasing modes. Another 
quantity which is useful for their description is the fractional comoving energy 
perturbation in each scalar field component 

Sef = m 3 + + 3ff^ - = d_ (fyA _ 

(c) 

where 5sj coincides with Bardeen's gauge-invariant quantity e m times the 
background energy density in the case of one scalar field (see [)I7|) and 
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relation (7) is used. Note also the following consequence of Eqs. (7-9) which 
is actually the Newton-Poisson equation in the cosmological case: 



N 

k 2 <$> = -AnGa 2 £ 8ef . (21) 

3=1 

In the long-wave limit k — > 0, the substitution of expressions (10,11) into 
(20) gives 0. This means that the small-A; expansion of Aj contains an ad- 
ditional k 2 multiplier in the case of both adiabatic modes, i.e., |Aj| <C |$| 
for them in this limit. On the other hand, |Aj| can be of the order of, and 
even much bigger than, |$| for isocurvature modes though the total comov- 
ing density perturbation J2f=i & £ f^ still contains the additional k 2 multiplier 
compared to $ as follows from (21). Substituting the expressions (16,17) 
valid during the multiple inflationary stage into (20) we get: 

Ai = 2diif + 87rG5^$> 3 = ^~,N. (22) 

3 

The next step is to determine the coefficients Ci, dj from amplitudes of 
quantum fluctuations of scalar fields generated during the inflationary stage. 
First, we invert (16,17) to obtain: 

(24) 

Further, using the above-mentioned possibility to add the same constant to 
all dj, we omit the last two terms in Eq. (24). Then all 5<pj in the r.h.s. 
of (23,24) have to be matched with quantum fluctuations of the scalar fields 
generated during the inflationary stage (we remind that this is a genuine 
quantum-gravitational effect). For all scalar fields being in the slow-rolling 
regime, l^f e //l = Wj \ ^ -^ 2 ' therefore, all mass- and $- dependent terms in 
Eqs. (8-9) may be neglected for k > aH, and even in the region k < aH but 
(t — tk)H(tk) H 2 (tk)/\H{tk)\ where t k is the time when a mode k crosses 
the Hubble radius during the inflationary stage, i.e., k = a{tk)H{tk) (it is 
in the latter region where the exact matching is performed). Then the Scft/s 
behave like massless uncoupled scalar fields in the de Sitter background. The 
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standard quantization gives the well-known result (see e.g. |Tj|): for k <C aH, 
the Fourier components of the fields are time independent ("frozen" ) and may 
be represented in the following form: 

<%(k) = ^Me,(k) (25) 
V 2k 6 

where ej(k) are classical stochastic Gaussian quantities with vanishing av- 
erage values < e,(k) >= and the correlation matrix < ej(k)e*,(k') >= 
5jj'8^ 3 '(h — k'). Note however that in the case of multiple inflation there may 



exist small effects |19j for which the approximation (25) is not sufficient and 
one has to take into account a small quantum correction to it reflecting the 
fact that the generated fluctuations are in a squeezed pure quantum state 
with a large but finite squeezing parameter r (the limit r — > oo is completely 
equivalent to (25)). As for scalar fields with large effective masses which 
are not in a slow-rolling regime, their fluctuations are negligible (apart from 
the case when they experience a non-equilibrium first-order phase transition 
during inflation which we do not consider here). 

By substituting (25) into (23,24), we get finally (we denote by C\ {k) the 
power spectrum of the stochastic quantity Ci(k) and use a similar notation 
for all stochastic variables, (/(k)/*(k')> = f 2 (k)5^(k - k')): 

<*« = -^fE^,, CfW-^E^, (26) 




k3 i V i 



9H 
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*( k ) = " ^ r , e > > ^) = 8ifc3V^' ^ = ^^ N ( 2? ) 



where all the time-dependent quantities in the r.h.s. are taken at t — tk- The 
result for Cf(k) coincides with that previously obtained in In the case 
of two scalar fields (j = 1,2), we reproduce the results of ref. where the 
notation C3 = d\ — d<i was used. The expressions (16,17,26,27) are the main 
results of this section. 



3 Matching to a post-inflationary era 

As was mentioned in the introduction, in the case of isocurvature modes 
we don't have general expressions like (10,11) for adiabatic modes, hence the 
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post-inflationary behaviour of isocurvature perturbations is not universal and 
depends on additional assumptions. In particular, there could be no such 
perturbations at all soon after the end of an inflationary stage. Therefore 
we will consider further a number of specific models in which they may be 
present even nowadays. The most natural way to achieve it is to assume 
that one of the inflaton scalar fields remains uncoupled from usual matter 
(baryons, photons, etc.) all the time since the end of the multiple inflationary 
stage up to the present moment, and that its particles or products of their 
decay (still uncoupled from usual matter) constitute today a part of the cold 
dark matter with a dust-like equation of state (p -C e). 

The non-decreasing mode of isocurvature fluctuations in a system of two 
uncoupled components consisting of dust-like matter on one hand, and ra- 
diation coupled to baryons on the other hand, can be characterized by a 
fractional comoving energy density perturbation in the dust-like component 

5 m = ^ = 5 i £-, Si = —1— « Aj , (28) 

that remains constant during the radiation-dominated era (see, e.g., |[20|-|[21| 
for reviews). Here Qi is the present-day density (in terms of the critical 
one) of that part of cold dark matter which is the relic from the inflationary 
era while Q m refers to all the cold dark matter (fij < fl m ). Note that we 
have for the number density perturbation, ^ — Si where n, is the number 
density of relics. VL tot = Q m + Q bar = 1 with great accuracy for cosmological 
models having an inflationary stage (the energy density of the cosmological 
term, if non-zero, should be added to fltot, too). After the transition to the 
matter-dominated stage at redshifts z ~ 10 4 , this mode produces a growing 
adiabatic mode of fluctuations which evolves, as usually, oc a(t) afterwards. 

Therefore, we have to relate Si at the radiation-dominated stage with Aj 
at the inflationary stage, as given in Eq. (22). We further specialize to the 
case of two inflaton scalar fields and replace the subscripts 1, resp. 2 by h 
(heavy), resp. I (light). Then Eqs. (16,17,22) take the following form which 
generalizes the results of ref. 0: 

$m - C (h) H I CM VlVI? ~ VhVl ' 2 (29) 

m - - Cl{k) lP + — (v h + v<)* ' (29) 

fit) = ^f + 2C 3(k )^V, (30) 
Oh ti Vh+ Vi 
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Ci(k) 



2C 3 (k) 



(31) 



H 



v h + vr 



A h (f) 



C7 s (k) \// 2 



(32) 



Ai(f) 



3 \4 + v; ' 
c 3 (k) 



(33) 



3 V h + Vi 



where C\ and C3 = d^ — di are given in Eqs. (26,27). Two essentially different 
cases may take place which we call the cases of heavy relics and light relics. 

1. Heavy relics. 

This case arises when the inffaton field that remains uncoupled from usual 
matter after inflation has an effective mass larger than H at the end of 
inflation. Then this "heavy" scalar field <ph is in the slow-rolling regime in 
the first part of inflation, but it goes out of this regime when H becomes less 
than the effective mass during inflation. Somewhat earlier, its energy density 
becomes much smaller than the total one. Let us take the simplest case 
where the effective mass is constant, so that the potential is Vh = m|0|/2, 
and G(j)\ ^> 1 at the early stages of inflation (for the field to be in the 
slow-rolling regime initially). Note for completeness that it is not possible to 
realize such a scenario for a steeper power-law potential 14, in particular for 
Vh = A0£/4, because then the effective mass remains smaller than H till the 
end of inflation. 

If Eh <C Stat, then irrespective of the fact whether the field <ph is in the 
slow-rolling regime or not, the right-hand side of Eq.(9) may be neglected for 
isocurvature modes. Now we need to solve this equation in the limit k 2 = 0. 
Note that then the left-hand sides of Eqs. (9) and (4) coincide in the case 
of a massive scalar field without self-interaction. So, one of the solutions is 
S(f>h oc 4>h{t) where <j>h(t) is the exact solution of Eq.(4) in a background driven 
by the other scalar field through Eq.(3). The other linearly independent 
solution can be found from the Wronskian condition but we don't need it 
here because, if (f>h is still in the slow-rolling regime, Eq.(30) may be applied 
which reads 5(j)h = 2C^H(j) h = —2C 3 m 2 h (f) h /3 for Vh <^ V. Now we use the 
constancy of the quantity 5<j>h(t)/(f)h(t) during a transition from inflation to 
the radiation-dominated stage. For rriht ^> 1 at the latter stage, the heavy 
field is in the WKB regime of oscillations with the frequency m h (see, e.g., 
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M for exact expressions). Averaging over the oscillations, we obtain 



Sh = 2 ^ = ~mlC 3 (k) . (34) 
9h 3 



Note that though Eq. (34) looks like Eq. (19.18) in j2(j], it is not exactly 
the same because it refers to a different quantity (a comoving energy per- 
turbation vs. an energy perturbation in the longitudinal gauge (6)), and we 
apply it in a different regime (at the radiation-dominated vs. the inflationary 
stage). From (27), the amplitude of the fractional energy and number den- 
sity perturbation in the relic dust component during the radiation-dominated 
stage follows: 

k 3 5 2 h (k) = 2H\t k ) (-L + ^) (35) 




where the first term inside the brackets in the r.h.s. of the last equation 
should be omittted if H(tk) < w,h- The power spectrum has a slope n « —3 
similar to that of $ 2 (fc) for adiabatic perturbations. 

To make a quantitative comparison between contributions of isocurvature 
and adiabatic modes to effects observable today, one should take into account 
that, due to the properties of the transfer function for isocurvature fluctua- 
tions in the CDM+radiation model, an isocurvature density fluctuation 5 m at 
the radiation-dominated stage produces the same adiabatic mode after tran- 
sition to the matter-dominated stage (a(t) oc t 2 ^ 3 ) as the initial adiabatic 
mode with $ = S m /5 for scales exceeding by far the present comoving scale 
corresponding to the cosmological horizon at the moment of matter-radiation 
equality R eq ~ 30h~ 1 Mpc, H Q = h ■ 100 km s^Mpc' 1 (Qbar is assumed to 
be small, too). For kR eq > 1, the equivalent amplitude of $ is even less. 
On the other hand, isocurvature fluctuations produce 6 times larger angular 
temperature fluctuations AT/T in the CMB at angles 9 > 30' for the same 
amplitude of long-wave density perturbations at the matter-dominated stage, 



i.e., AT/T = 25 m /5 vs. AT/T = $/3 for adiabatic perturbations p2 |, [^3 

(see also JZIJ-^TJ for reviews). Due to the latter reason, it has been 
long known that isocurvature fluctuations with a flat (n = —3) initial spec- 
trum cannot be responsible for the observed large-scale structure and AT/T 
fluctuations in the Universe. 

For a power-law VJ with the last part of inflation driven by the light scalar 
field, VI > V{ in the region Vh ~ VJ where the transition from heavy to light 
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scalar field domination of the total energy takes place. Then the second 
term inside the brackets in (35) is the dominant one, while still much smaller 
than $> 2 (k) = 9Cf(fc)/25. Therefore isocurvature fluctuations, if present, are 
less than adiabadic ones in double inflationary models in the region around 
the break in the perturbation spectrum due to a transition between the two 
phases of inflation, and their possible presence changes nothing regarding the 
confrontation of these models with observational data (see, e.g. Ul2| ). But 
on much smaller scales, when the first term inside the brackets in (35) is 
dominant, isocurvature perturbations become much larger. In that case 

kX(k) = 2J f(h) ■ (36) 

Alternatively, this result may be obtained very simply by considering the 
heavy field as a test field in the de Sitter background and using the expres- 
sions (25,34). The spectrum (36) grows with k (because 4>h quickly decreases 
with t) until the point H{tk) ~ rrih is reached, after that it falls abruptly. 

If, e.g., Vi = \mf(f)f and mi <C m^, then, using Eqs. (2.10-2.15) of ref. |§ 
and Eq. (35), we get 

kX(k) ^(_^_y < + ^ yi k>>h (37) 

where kb is the location of the break and so ^> 1 is the number of e-folds 
during the second phase of inflation driven by the light scalar field (so ~ 60 
to account for observational data, see Jl2|). The expression (37) is derived 
under the approximation m\/mf < sq which corresponds to the absence of 
a power-law intermediate stage between the two phases of inflation (double 
inflation without break, according to the terminology of |J), a more suitable 
condition perhaps is the absence of oscillations or a smooth transition in 
the spectrum, which will be the case for m^/m/ < 15 2l| or ml/mf < 4sq. 



In the opposite case of double inflation with a break, there is no growth of 
isocurvature fluctuations at small scales. Note that the effect of growth in the 
isocurvature perturbation spectrum was previously noticed from numerical 
calculations for a similar model in 0. 

For not too small scales when ln(k/kb) Sq, the first term in the spec- 
trum (37) is power-law like with a small exponent: k 3 5l(k) oc (k/kb) a , a = 
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m\/mfso. The perturbations reach their maximum on short scales for which 
s — ha.(k/kb) ~ ml/mf, due to the disappearance of the first term in (37) 
for H{t k ) < m h , its value being 



. (38, 

It is interesting that for VGmi ~ 10~ 6 , s ~ 60 the maximal value of the 
quantity (38) as a function of rrih/mi, though still smaller than unity, is not 
far from it (it is reached for rrih/mi ~ 5). Hence, such a model can be used 
to produce a significant number of primordial black holes with rather small 
masses (for a review of observational upper limits on the number density 
of PBHs see, e.g., pB"|-[p7|). Then, however, it cannot explain the observed 



large-scale structure and AT/T fluctuations in the Universe because m^jmi 
is required to be ~ 12 — 14 (and certainly more than 8) for this aim, see 



2. Light relics 

In this case, vr? e ^ <C H 2 for one of the inflaton scalar fields during the whole 
inflation. To avoid this "light" field to be dominating during the last part of 
inflation and after its end, we have to assume that its energy density ei <C e to t 
during inflation, too. Then we have from Eq. (31): 

8(j>i = -2C 3 H^ = ho.Vl . (39) 

If C3 from Eq. (27) is substituted into this expression and it is assumed that 
V{ <C Vfr, the standard expression (25) for fluctuations of a test scalar field on 
a de Sitter background arises once more. Let Vi = ^mfcj) 2 , then we may repeat 
the derivation made in the previous subsection to get the expression (36) 
(with the index "h" changed to "1") for the fractional density perturbation 
in the light relic component. Now <pi is practically constant during inflation 
(and less than G -1 / 2 to avoid a second inflationary phase), so the spectrum 
is falling with k. Isocurvature perturbations may be larger than adiabatic 
ones if <pi is small enough, but this does not lead to interesting cosmological 
models for a smooth potential Vh satisfying the slow-rolling conditions due 
to the reason mentioned in the previous subsection in connection with the 
n = —3 initial isocurvature perturbation spectrum. 
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3. Intermediate relics 

Let us briefly consider the case of an inflaton field which remains uncoupled 
from usual matter after inflation, dominates during the first phase of inflation 
(so we call it "heavy") and has the quartic potential Vh = A/ l 0^/4. During 
the last part of inflation, s < s , where s is the number of e-folds measured 
from the end of inflation and Sq ^> 1 is the moment when Vh = V, <C e tot , 
however m\ e ff = 3\ h (f) 2 h <C H 2 (and less than the effective mass of the other 
scalar field, too). So, this initially "heavy" inflaton becomes "light" in the 
last part of inflation. That is why we call this case the intermediate one. 

Then, for s <C so, S(ph = ^C^Htph = —2C^\h4>h{t)/?> as in the case of mas- 
sive relics. Here the quantity 5(fih/(fih is not constant during the last period 
of inflation and the transition to the radiation-dominated stage. Therefore, 
an exact matching (as it was done in the subsection 1) is not possible, but 
we may make a matching by order of magnitude using the fact that Sh at the 
radiation-dominated stage is of the order of 54>h/4>h at the end of inflation. 
Using (30,34), we get: S^(k) = const ■ C%(k)m\ e j j(t /) where tf is the moment 
when inflation ends and const = 0(1). 

a) V = m 2 (fi 2 /2 , \ h > Gmf . 
Then <p\ = mf/X^ ln(s /s) during the last period of inflation. Therefore, 
m h,eff(tf) = 3m^/lns , and we arrive at the following result 



where s(k) = s(t fc ) = \n(kf/k) 3> 1 and kf = a(tf)H(tf). The spectrum 
is approximately flat. Though it has a smooth maximum at s = So e ~ 3 ) 
this maximum is not as strongly pronounced as in the case of massive heavy 
relics. The amplitude grows oc {k/k^) 1 ^ for s Q — s <^ s where kb is the inverse 
comoving scale corresponding to the first horizon crossing at the moment s . 




(40) 



where const 



(9(1). In particular, for s <^ s , 




(41) 
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b) V x = A^f/4 , \ h > A, . 
Now <f>1 = \i(pj/2\h ln(s /s) during the second phase of inflation. Thus, 
m h,eff(tf) ~ A^/Glnso and 

*■«<*) = const • tfft) ( + ^ (42) 

where const = 0(1). In particular, for s <C s , 

kV d k) = const (43) 

The spectrum is approximately flat and grows slightly towards large /c's. 
Once more, its amplitude grows oc (fc/fcj,) 1-5 for s — s ^ s o- 
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